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Abstract. The most massive black holes, lurking at the centers of large galaxies, 
must have formed less than a billion years after the big bang, as they are visible to- 
day in the form of bright quasars at redshift larger than six. Their early appearance 
is mysterious, because the radiation pressure, generated by infalling ionized mat- 
ter, inhibits the rapid growth of these black holes from stellar-mass black holes. It is 
shown that the supermassive black holes may form timeously through the accretion of 
predominantly degenerate sterile neutrino dark matter onto stellar-mass black holes. 
Our symbiotic scenario relies on the formation of, first, supermassive degenerate ster- 
ile neutrino balls through gravitational cooling and, then, stellar-mass black holes 
through supernova explosions of massive stars at the center of the neutrino balls. 
The observed lower and upper limits of the supermassive black holes are explained 
by the corresponding mass limits of the preformed neutrino balls. 



Keywords: dark matter, neutrino properties, massive black holes 
ArXiv ePrint: ??? 



Author to whom any correspondence should be addressed 



A symbiotic scenario for the rapid formation of supermassive black holes 



2 



1. Introduction 

Supermassive black holes of ~ 3 x 10 9 M© p] are already present at redshift z = 6.42 or 850 Myr after the big 
bang, as is evident from the recently discovered [2] quasar SDSS J114816. 64+525150. 3. Unlike the galaxies, 
which are thought to assemble through hierarchical mergers of smaller galaxies [3] , the supermassive black 
holes form anti-hierarchically, i.e. the larger black holes prior to the smaller ones. This is supported by the 
observation [4] that the number of brighter quasars peaks at z ~ 2, while that of lower-luminosity active 
galactic nuclei peaks at z ~ 0.6. There is only a narrow window of about 485 Myr between the appearance 
of the first stars at z <~ 11, or 365 Myr after the big bang when reionization of the universe begins [5], and 
the appearance of the first quasars [5] at z — 6.42 when reionization ends, during which at least some of the 
most massive black holes must have formed. 

It has been known for more than three decades [6] that the growth of a black hole from the lOM© to the 
3xlO 8 M0 scale, through the accretion of baryonic matter radiating at the Eddington limit, with an efficiency 
between e ~ 0.1 and e ~ 0.2, would take at least 1.6 Gyr. This seems to preclude the existence of a ten 
times larger black hole at z = 6.42 that must have been produced in only ~ 0.5 Gyr. Such a short accretion 
time can, in the absence of frequent black hole mergers, only be achieved by drastically increasing either 
the ambient baryonic matter density, or the mass of the seed black hole, or both [5]. The baryonic matter 
density would have to be much larger than the average matter density within 1 pc of the Galactic center 
of ~ lO 7 M /pc 3 , which is already uncomfortably large. Alternatively, increasing the seed black hole mass 
by several orders of magnitude is not attractive either, because, in spite of the ongoing intensive search, the 
evidence for intermediate mass black holes, with masses between 20 M Q and 10 6 M Q , is rather weak and 
controversial [7] . Regardless of whether the supermassive black holes were produced through the accretion 
of baryonic matter onto stellar-mass black holes, or through binary mergers starting from stellar-mass black 
holes, these would inevitably have left a trail of observable intermediate mass black holes, which ran out of 
baryonic matter or black hole supplies sometime and somewhere in this universe. 

Thus apart from the easily detectable upper limit of M max ~ 3 x 1O 9 M , there also seems to be a lower 
mass limit of M m \ n ~ 10 6 M Q of the supermassive black holes. In this context, it is important to note that 
the mass of the black hole at the center of M87 [B] is equal to that of the earliest quasar [5] observed at z 
= 6.42. As these are both archetypical examples of the most massive black holes at vastly different epochs 
of the universe, one may infer that M max has not increased significantly from z = 6.42 to z = 0, or from 
0.85 Gyr to 13.7 Gyr after the big bang. Indeed, although Eddington-limited baryonic matter accretion is 
essential for seeing the quasars, this is a transient phenomenon with an estimated total lifetime of a few tens 
of Myr [3] , that presumably contributes very little to the final mass of the most massive black holes today. 

In summary, it seems that a consistent theory of the formation of the supermassive black holes should be 
able to explain both observational facts, the nearly time-independent upper and lower mass limits of the 
supermassive black holes, as well as the early formation of the most massive black holes around z = 6.42. 
As both the baryonic matter accretion, as well as the black hole merger scenarios, do not seem to provide us 
with a satisfactory description of these observational facts, we are led to ask the pertinent question, whether 
an alternative scenario, based on the accretion of mainly dark matter, may do better. Galaxies are, indeed, 
dominated by dark matter, and as part of this dark matter is concentrated in the galactic centers, it may very 
well have contributed to the formation of the supermassive black holes. Of course, in order to make definite 
predictions, one needs to focus on a well-defined and consistent dark matter candidate. Thus, in section 2, 
we discuss the physical and cosmological properties of our sterile neutrino dark matter candidate, while in 
section 3 we explore an astrophysical consequence of this dark matter scenario, namely the formation of sclf- 
gravitating supermassive degenerate sterile neutrino balls, with masses between 1O 6 M0 and 3 x 10 9 Mq. We 
then discuss how these preformed neutrino balls convert efficiently and anti-hicrarchically into supermassive 
black holes, using stellar-mass seed black holes as catalysts. In section 4, we describe the dynamics of the 
accretion of degenerate sterile neutrinos onto a black hole in a simple nonrelativistic Thomas-Fermi field 
theory, based on the Lane-Emden equation, while our conclusions are presented in section 5. 
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2. Physics and cosmology of sterile neutrinos 

Recently, a set of three right-handed sterile neutrinos has been consistently embedded in a renormalizable 
extension of the minimal standard model of particle physics, dubbed the vMSM [9]. While two of these 
sterile neutrinos are unstable, having masses in the 1 GeV/c 2 to 20 GeV/c 2 range, the third one, with mass 
around 10 keV/c 2 , is a promising quasi-stable dark matter candidate. This light sterile neutrino interacts 
with the standard model particles only through its tiny mixing with the active neutrinos. The bulk part of 
the light sterile neutrinos is produced ~ 2.3 [is after the big bang, at temperatures 

well ahead of the quark-gluon and chiral restoration phase transitions, through incoherent resonant [10] and 
non-resonant [11) scattering of the active neutrinos. For a wide range of the parameters of the ^MSM, the 
sterile neutrinos are generated out of thermal equilibrium, yielding a sterile neutrino mass fraction of the 
total mass-energy of this universe which is consistent with that of nonbaryonic dark matter. For instance, 
with an initial lepton asymmetry 

= n^n^ = 10 _ 2 (2) 
n 7 

a mixing angle of the sterile neutrino v s with the active neutrino v e given by sin 2 $ = 10~ 13 , and a sterile 
neutrino mass m = 15 keV/c 2 , the mass fraction of the sterile neutrinos produced in the early universe is 
indeed Q„ a ~ 0.24 10J, i.e. equal to that of the nonbaryonic dark matter, derived from the WMAP data [5]. 

The same tiny mixing angle, which prevents the thermal equilibration and thus the overproduction of the 
sterile neutrino in the early universe, also renders the sterile neutrino quasi-stable [12] , with a lifetime much 
larger than the age of the universe. In fact, with our choice of the parameters, sin 2 $ — 10~ 13 and m = 15 
keV/c 2 , the sterile neutrino decays with a lifetime of 

T {y B -> v e vv) = - = 1.21 x 10 19 yr , (3) 

Gpm b sin v 

predominantly into a v e and a neutrino-antineutrino pair, and iJi, carrying the flavours i = e, /i or t [12]. 
These standard neutrinos are all virtually unobservable because their energy is too small. But there is also 
a subdominant radiative decay mode, with a branching ratio 

; = — = 0.784 x 10" 2 (4) 



into a potentially observable photon and a v e [12j . However, due to the smallness of its partial decay width 
of 

[t(v s -> VelT 1 = 0-649 x 10" 21 yr _1 , (5) 

these photons of energy mc 2 /2, which may be the "smoking gun" of the sterile neutrino, are difficult to 
observe as well. In fact, for the chosen model parameters, a sterile neutrino dark matter concentration of 
mass M has a luminosity of merely 

Mc 2 

Lx = T-7 \ = 1-84 x 10 25 (M/M ) erg/s , (6) 

in photons of mc 2 /2 = 7.5 keV energy. Thus the best places to look for these photons are the diffuse 
extragalactic X-ray background, as well as the X-rays emitted by large galaxy clusters, low-surface-brightness 
and dwarf galaxies that are dominated by nonbaryonic dark matter [13] . For an initial lepton asymmetry 
L V(S ~ 10~ 10 , which is of the order of the baryon asymmetry [5] 

S = I^6xlO~ 10 , (7) 

n 7 
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mainly non-resonant neutrino scattering contributes to the production of sterile neutrino dark matter. These 
sterile neutrinos inherit a nearly thermal energy spectrum from the active neutrinos [11] , which allows them 
to play the role of warm dark matter in the large-scale structure of the universe, the clusters of galaxies and 
the galactic halos. They may also erase the undesirable excessive substructure on the galactic scales. 

The initial lepton asymmetry does not need to be of the same order of magnitude as the baryon asymmetry. 
However, for a larger initial lepton asymmetry, like L Ve ~ 10 -2 , there is, in addition to non-resonant 
neutrino scattering, also resonant or matter-enhanced neutrino scattering contributing to the production of 
dark matter [TUJ. The latter yields cool sterile neutrinos that have a distorted quasi-degenerate spectrum, 
with an average energy of about two-thirds of that of the warm sterile neutrinos, due to the resonant 
Mikheyev-Smirnov-Wolfenstein (MSW) oscillations [14] . 

A relatively large initial lepton asymmetry of L Uc ~ 1CP 1 to 1CP 2 is also what seems to be required, to bring 
the observed light element abundances in line with the number of three active neutrinos at nucleosynthesis 
|15j . Sterile neutrinos may as well be responsible for the pulsar kicks of up to ~ 1600 km/s, which magnetars 
acquire in supernova explosions [16j . Thus our sterile neutrino meets all the constraints which an acceptable 
dark matter particle must fulfil [10j.|13j. but it has several remarkable additional properties that make it a 
rather unique candidate for dark matter. 

3. A symbiotic black hole formation scenario 

For our model parameters m — 15 keV/c 2 , sin 2 i? = 10~ 13 and L Ue = 10~ 2 , cool (or resonant) dark matter 
dominates over warm (or non-resonant) dark matter by a factor of about three [IS] . The cool sterile neutri- 
nos become nonrelativistic ~ 22 min after the big bang, well after nucleosynthesis, and they begin, together 
with the warm sterile neutrinos and baryonic matter, to dominate the expansion of the universe ~ 79 kyr 
after the big bang, well ahead of recombination. Thus the primordial density fluctuations of sterile neutrino 
dark matter have enough time to grow nonlinear and form degenerate sterile neutrino balls [17] . through a 
process called gravitational cooling [18 , prior to the appearance of the first quasars, 850 Myr after the big 
bang. This collapse process may start ahead of reionization, perhaps as early as ~ 320 Myr after the big 
bang. Initially, the free-falling sterile neutrino dark matter, dominating baryonic matter by about a factor of 
six [5], drags the baryonic matter along towards the center of the collapse. The baryonic gas will get heated, 
reionized and evaporated, but the free fall of the neutrinos will not be inhibited by the Eddington radiation 
limit. Eventually, the quasi-degenerate sterile neutrino dark matter hits, ~ 640 Myr after the big bang, the 
degeneracy pressure, bouncing off a number of times, while ejecting a fraction of the dark matter at every 
bounce. The neutrino ball finally settles in a condensate of degenerate sterile neutrino matter at the center 
of the collapsed object, as has been shown in calculations based on time-dependent Thomas-Fermi mean 
field theory [IS]. 

The smallest mass that may collapse is the mass contained within the free-streaming length at matter- 
radiation equality, ~ 79 kyr after the big bang. For m = 15 keV/c 2 , this free-streaming mass is M warm ~ 
7 x 10 6 M Q in the case of warm (or non-resonant) sterile neutrinos, and M coo i ~ 2 x 1O 6 M in the case of 
the dominant cool (or resonant) sterile neutrinos 10J . As part of the neutrino dark matter is ejected during 
the collapse process, the minimal mass of a degenerate sterile neutrino ball may be somewhat smaller than 
■Wcool) perhaps M m j n ~ 10 6 M Q , consistent with the lower mass limit of the observed supermassive black 
holes. 

The maximal mass that a self-gravitating degenerate neutrino ball can support gravitationally, is the 
Oppenheimer-Volkoff (OV) limit Qjj] 



For Tji — 15 keV/c 2 and spin degeneracy factor g = 2, this mass M max = 2.789 x 10 9 M Q is consistent with 
that of the most massive black holes observed in our universe QQ,P]. As such a supermassive degenerate 
neutrino ball has a radius of only 4.45 Schwarzschild radii ;19. , it is almost a black hole. Thus the maximal 
mass scale of these objects may be linked to the existence of a sterile neutrino of ~ 15 keV/c 2 mass, in a 
similar fashion as the maximal mass scale of the neutron stars is linked to the effective mass of the neutron 



M, 



max — 



0.5430 Mp 




(8) 
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Since the gravitational potential in the interior of a neutrino ball is nearly harmonic, these objects, in 
particular those near the upper mass limit, are ideal breeding grounds for stars of mass M > 25M . As 
soon as such a central massive star is formed from a collapsing molecular hydrogen cloud that was attracted 
to the neutrino ball, it may be kicked out through close binary encounters with intruding stars. However, 
before that happens, a star of 25 Mq will evaporate large portions of its hydrogen and helium envelope 
and become a Wolf-Rayet star. And about 3 Myr after its formation, the star undergoes a core collapse 
supernova explosion of type Ic, leaving a black hole of 3 to 4 Mq at the center of the neutrino ball. Some 
of these most massive supernova explosions, occurring in high-mass neutrino balls, may be observable as 
long-duration 7-ray bursts [21]. As in contrast to pulsars, black holes do presumably not acquire "black hole 
kicks" during a supernova explosion, the velocity of the stellar-mass black hole will be small compared to the 
escape velocity from the center of the neutrino ball. For a ball of 3xlO 6 M0 mass and 25 light-days radius 
[T7] . consisting of degenerate sterile neutrinos of 15 keV/c 2 mass and degeneracy factor g — 2, the escape 
velocity from the center is 1700 km/s, while for a ball of the same sterile neutrinos at the OV-limit, with 
2.8xl0 9 Mq mass and 1.4 light-days radius, the escape velocity from the center is the velocity of light [19|. 

The supernova explosion of the massive star, giving birth to a stellar-mass black hole at the center, sparks the 
rapid growth of the black hole through nearly radiationless, and therefore, Eddington-unlimited accretion 
of mainly degenerate sterile neutrino dark matter from the surrounding neutrino ball, until the supplies 
dry up. In this symbiotic scenario, the anti-hierarchical formation of the bright quasi-stellar objects and 
low-luminosity active galactic nuclei may be explained by the fact that the escape velocity from the center 
of a 3 xlO 6 M0 neutrino ball, is 176 times smaller than that of a 3 xl0 9 Ai© neutrino ball. The low-mass 
neutrino ball may, therefore, have difficulty capturing a molecular hydrogen cloud that is able to produce 
a massive star. In particular, a low-mass neutrino ball may experience a large number of unsuccessful 
attempts, leading to ordinary low-mass stars, or neutron stars after supernova explosion, prior to delivering 
the expected stellar-mass black hole. These unwanted stellar-mass objects will eventually be ejected from the 
neutrino ball through pulsar kicks or close binary encounters with intruder stars from the surrounding star 
cluster, thus clearing the scene for the next attempt at forming this stellar-mass black hole. The randomness 
of this process may very well delay the formation of a stellar-mass black hole at the center of a low-mass 
neutrino ball by several Gyr, while a neutrino ball at the top of the mass scale may easily deliver the stellar- 
mass black hole on its first attempt in less than 10 Myr. Of course, these sketchy ideas will have to be tested 
in realistic numerical simulations. However, if this scenario is correct, some low-mass neutrino balls may 
still be around at some galactic centers. For instance, a 10 6 Mq neutrino ball would reveal itself through its 
X-ray emission of 2 x 10 31 erg/s at 7.5 keV, for our model parameters. 



4. Accretion of a neutrino halo onto a black hole 

A particle that is initially at rest at the surface of a 3 x 10 6 M Q neutrino ball reaches the center in the 
free-fall time Tp ~ 35 yr. This is also the time frame in which the accretion process onto the central black 
hole reaches a steady-state flow. In the steady-state approximation, the flow is governed by Bernoulli's eq. 

0(r) + ^(u 2 {r) + v 2 F {r)) = <j){r H ) = const . (9) 

Here u(r) is the flow velocity of the infalling degenerate sterile neutrino fluid, vp(r) its Fermi velocity, <p(r) 
the gravitational potential and r# the radius of the halo. Assuming that the flow makes the gravitational 
potential </> extremal for all values of the radius r, with respect to variations that satisfy the constraint of 
mass conservation, 

m 4 g v% u ,„ „. 

pu = 2f * 3 = const , (10) 

we obtain 

u 2 (r) = ^v%(r) = 4(r) , (11) 
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which means that the inflow of the sterile neutrinos is trans-sonic, i.e. it flows at the local velocity of sound 
cs(r). Thus Bernoulli's eq. is now 



y F (r) = <f>(r H ) - 0(r) = GM Q 



v(x) 
bx 



(12) 



which defines the quantity v(x). Using Poisson's eq., one can readily verify that v(x) fulfils the Lane-Emden 
eq. 



1 d 2 v(x) 
x dx 2 



v(x) 



3/2 



provided the length scale is 



b = 



3nh 3 



2/3 



(13) 



(14) 



3 V4V2<?to 4 G 3 / 2 a4 /2 J 

Thus for m = 15 keV/c 2 and g = 2 we have b = 2.587 lyr. The total mass enclosed within a radius r = bx 
is [22] 

M(x) = M [v(x) - xv '(x)] . (15) 

Various solutions v(x) of the Lane-Emden eq.(13), all having the total mass M = 2.714 M©, are shown in 
Fig.l. There are three distinct classes of solutions. The M-solutions exhibit a central point mass Mp = 



v(x) 

o 

2.5 :: : 



1.5 v 



0.5 - 




1 2 



3 4 5 6 



Figure 1: Various solutions of the Lane-Emden equation, all having total mass M = Mp+Mn — 2.714 Mq. 
The solid, dashed and dash-dotted lines represent the E— , F— and M— solutions, respectively. 



Mqv(0), surrounded by a self-gravitating degenerate sterile neutrino halo. The F-solutions describe shells of 
self-gravitating degenerate neutrino matter that are gravitationally unstable. The E-solution, with v(0) = 
and v'(0) = 1, stands for a pure neutrino ball with Mp = Mqw(O) = 0. Our focus is on the M-solutions of 
the Lane-Emden eq.(13) because these describe, in the steady-state approximation, the various stages of the 
accretion history of a black hole surrounded by a degenerate sterile neutrino halo. For instance, decreasing 
the halo radius ru — bxn, while keeping the total mass M fixed, causes the central point mass Mp to 
increase, as seen in Fig.l. The solutions of the Lane-Emden eq.(13) with Mp = M@v(0) > 0, for arbitrary 

mass M can be obtained noting that, if v(x) is a solution, v (x) — A 3 v(Ax) with A > is a solution as 
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well. Thus all the masses and radii scale as M — A 3 M and r = r/A [22]. The mass accretion rate into a 
sphere, containing the mass Mq within a radius tq from the center, is thus given by 

— £ = Airr c p(r c )u(r c ) (16) 

or 

dM c y/Zgwt&Ml 2 f(p) 

~w = — 2^ — Hxc)] = Mc ' (17) 

where we have introduced the universal time-scale 

2irh 3 7 2 Zl5keV\ 4 . , 

= 1.488 xlO 7 - — yr . (18) 



V3gm 4 G 2 M Q g \ mc 2 

The shut-off parameter 

M = ^ 2 (19) 

[v{x c ) - xcv'{xc)] 

is a function of the mass ratio \i = Mc/M. We now choose rc — bxc to be the radius at which the escape 
velocity reaches the velocity of light, i.e. xc is given by 

J- = ^-«'(*h) , (20) 

Xq X C 

where r Q = bxQ is the Schwarzschild radius of the sun with 



2GM a „ /o\2/3 

u — 1 1(17 \s m-u / a 



2 \ 8/3 



* e = -j^. 1.207X10-"^ . (21) 

As for M <C M max and \i < 1, we may approximate w(icc) ~ u(0) and /(/i) ~ 1, eq.(17) agrees well with 




Figure 2: Various shut-off parameters /(/i) as a function of the mass ratio [i = Mc/M . The neutrino ball 
masses vary from M = 10 6 M Q for the box function, through M = 10 7 , 10 8 to 10 9 M for the curve with the 
largest peak value. 
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standard Bondi accretion theory [23] 
hole is given by 



Integrating eq.(17) using this approximation, the growth of the black 
M c (0) 



M c (t) 



1 - t/T A 



(22) 



yielding an accretion time-scale of ta = tMq/Mc{0)- During the accretion process, both Mq and xc grow, 
while xh shrinks as a function of time, eventually causing xc and xh to converge and v(xc) to vanish. The 
shut-off parameter /(/i), shown in Fig. 2 as a function of the mass ratio (x = Mq/M, is for a neutrino ball of 
mass M = 10 6 Mq a simple Heavyside function. As M increases towards M max , this curve starts deviating 
from the simple box form, thus signalling the breakdown of our nonrclativistic theory. Since, for M -C M max , 
the black hole growth is approximately given by the Bondi formula, Mp oc Mq, we expect the mass growth 
curves to match the Bondi solution closely, with the attractive feature that the solution of eq.(17) eventually 
brings the growth to a halt. The mass ratio fi is shown in Fig. 3 as a function of time in units of the accretion 
time-scale ta- For low-mass central black holes, these curves are indistinguishable from the Bondi solution, 
represented by the dotted curve. Here we start with /i(0) = 0.1 at t = to illustrate the differences for halo 
masses between 10 6 and 10 9 M Q . The universal time-scale is r = 14.88 Myr, and the accretion time-scales 
for a 3x 10 6 M© neutrino ball onto 3 and 4 Mq seed black holes are, therefore, ta = 4.96 Myr and ta = 3.72 
Myr, respectively. 




Figure 3: The growth of a central black hole as a function of time, in units of the accretion time ta, for 
halos of between 10 6 and 10 9 M Q , as in Fig. 2. For fj, < 0.1, the growth curves match the Bondi growth curve 
too closely to distinguish. 



5. Conclusions 

In summary, the neutrino balls are almost entirely swallowed by the seed black holes in an accretion time 
scale ta < 5 Myr, thus converting these rapidly into supermassive black holes with negligible residual sterile 
neutrino halos. The most massive black holes may, therefore, form between 650 Myr and 840 Myr after the 
big bang. Although the nonrelativistic Thomas-Fermi mean field theory breaks down for M ~ M max and 
r H ~ r Ci we expect these basic results to persist in a relativistic theory of the accretion process. 

There are three main features which distinguish the accretion of degenerate sterile neutrino dark matter 
from that of baryonic matter, playing a decisive role in the rapid growth of a stellar-mass black hole to 
the supermassive scale. Firstly, in contrast to the clumping of baryonic matter, neither the neutrino ball 
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formation nor the neutrino halo accretion onto a black hole is inhibited by the Eddington radiation limit. 
Secondly, the matter densities of the degenerate sterile neutrino balls are much larger than those of any 
known form of baryonic matter having the same total mass, leading to much faster growth of the black holes 
through neutrino dark matter. Thirdly, the preformed degenerate sterile neutrino balls have, for m ~ 15 
keV/c 2 and g = 2, masses in the same range as the observed supermassive black holes, which sets a natural 
limit to the growth of the black holes. 

We, therefore, conclude that supermassive neutrino balls, with stellar-mass black holes at their center, indeed 
offer an intriguing symbiotic scenario, in which baryonic matter conspires with degenerate sterile neutrino 
dark matter, to form these galactic supermassive black holes, with masses between 10 6 M© and 3 x 10 9 M Q , 
rapidly and efficiently. 
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